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ABSTRACT 



In this Letter, we report the discovery of a new bright radio transient in M82. Using the Very Large Array, we observed the nuclear 
region of M82 at several epochs at 22 GHz and detected a new bright radio source in this galaxy's central region. We find a flux 
density for this flaring source that is ~ 300 times larger than the upper limits determined in previous observations. The flare must have 
started between 2007 October 29 and 2008 March 24. Over the past year, the flux density of this new source has decreased from ~100 
mJy to ~ 1 1 mjy. The lightcurve (based on only three data points) can be fitted better with an exponential decay than with a power 
law. Based on the current data we cannot identify the nature of this transient source. However, a new radio supernova seems to be the 
most natural explanation. With its flux density of more than 100 mJy, it is at least 1.5 times brighter than SN1993J in M81 at the peak 
of its lightcurve at 22 GHz. 
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1. Introduction 

M 82 is a nearby (3.6 Mp c based on a Cepheid distance to M81 
by iFreedman et al] 1 19941) irregular (10) galaxy with a very ac- 
tive starburst in its nuclear region. It harbors many bright su- 
pernova remnants in its ce ntral region, which have been stud- 
ied extensively for decades (Muxlow et al. 1994; Beswick et alj 
l2006; Fenechetal. 2008). van Bureii &_Greenhouse (1994) es- 
timate a supernova rate of 0. 1 yea r~'. However , a new radio su- 
pernova has not been discovered. ISinger et al.l (|2004|) reported 
a supernova in M 82 (SN2004am) that was classified as type-II 
jMattila et al. 2004). b ut it h as not been detected at radio wave- 
lengths ([Beswick et al.j 2004. 

iKronberg & Sramekl (Il985h and lKronberg et al.1 (l2000l) mon- 
itored the flux densities of 24 radio sources in M82 from 
1980 until 1992. Most sources (75%) remained surprisingly 
constant. There is some controversy about how the fluxes of 
these compact radio sources can be stable. Models of super- 
nova renmants_ex£andin^ dense medium may explain 
fliis jChevalier & F ransson" 200ll) . ISeaquist & StankovicI (l2007h 
argue that the radio emission could arise from wind-driven 
bubbles. Studying the evolution of a young source could be 
very important for understanding these models. The strongest 
source, 41.9-1-58.0, sho ws an exponential de cay with a decay 
rate of t^=11.9 years ( Kronberg.etaD |2000|). Another source, 
41.5-1-597, which was detected in 1981 with a flux density of 
~10 mJy, faded within a few months to a flux densi ty below 1 
mJy (IKronberg & Sramek|[T985t IKronberg et aI.ll200Ql) . The na- 
ture of this strongly variable source was never clarified. 

Radio supernova are rare events. So fa r only about two 
dozen have been detected dWeiler et al.ll2002l) and most of them 



were quite distant and rather weak. This makes it difficult to 
study them in great detail. One notable exception is SN1993J 



(ISchmi dt et all 1 1 9931) in M81, which has been studied exten- 
sively ( Marcaide etdIT 997. 2009; Bietenholz et al. 2001.|2003t 
iPerez-Torres et aI.ll200U i2002i; ,BarteI et al.ii2002t i2007i) . Thus, 
the detection of a new nearby supernova would be highly desir- 
able. 

M82 is part of the M8 1 group of galaxies and it s hows clear 
signs of tidal interaction with M81 and NGC3077 (lYun et al.l 
Il994l) . This makes the M81 group an ideal system for study- 
ing galaxy interaction in great detail. With this motivation, 
we have initiated a project to measure the proper motions of 
M81 and M82 with VLBI astrometry. We are observing M81*, 
the nuclear radio source in M81, bright water masers in M82, 
and three compact background quasars. Based on our experi- 
enc e with mea surements of proper motions in the Local Group 
(Brunthaler et al. 2005b, 2007) . we expect a detection of the tan- 
gential motions of M81 and M82 relative to the Milky Way 
within a few years. So far, we have observed at three epochs 
at 22 GHz with the High Sensitivity Array (including the Very 
Long Baseline AiTay, the Very Large AiTay, and the Greenbank 
and EfFelsberg telescopes). Here, we report the detection of a 
new transient source in M82 based on the data from the NRACQ 
Very Large Array (VLA). 



' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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2. Observations and data reduction 

M82 was observed with the Very Large Array (VLA) as part 
of the High Sensitivity Array observation under projects BB229 
and BB255 on 2007 January 28, 2008 May 03, and 2009 April 
08. The total observing time in each epoch was 12 hours. We 
used M81* as phase calibrator and switched between M81*, 
M82, and 3 extragalactic background quasars every 50 seconds 
in the cycle M81* - 0945+6924 - M81* - 0948+6848 - M81* 
- M82 - 1004+6936 - M81*, yielding an integration time of 
~100 minutes on M82. We observed with two frequency bands 
of 50 MHz, each in dual circular polarization. The data reduc- 
tion was performed with the Astronomical Image Processing 
System (AIPS) and involved standard steps. On 2007 January 
28 and 2009 April 08, we used 3 C48 as flux density calibrator 
M81* was used as a gain and phase calibrator and one round 
of phase self-calibration was performed on M82. Unfortunately, 
no flux density calibrator was observed on 2008 May 03. Here, 
we assumed a flux density of 150 mJy for the highly variable 
source M81*. This value was chosen si nce it was in the range 
of typical values at cm wavelenghts (e.g. lBrunthaler et"aLll2006l) 
and it yielded flux densities for 0945+6924, 1004+6936, and 
44.0+59.6 in M82 that were consistent with their flux densities 
at the other two epochs. 

M82 was also observed on 2008 March 24 with the VLA at 
22 GHz for 10 minutes (~6 minutes integration time on M82) in 
spectral line mode. 3C 48 was used as flux density calibrator, and 
1048+717 was used as gain and phase calibrator. A total band- 
width of 9. 18 MHz was observed. We also analyzed the archival 
VLA data of 2007 October 29 at 4.8 GHz, which is the latest 
available data before our observations. Here, 3C286 was used 
as flux density calibrator and 1048+717 was used as phase cali- 
brator. Eighteen minutes of integration time on M82 was spread 
over 1.5 hours. 



3. Results 

A new bright source was detected on 2008 May 03 (see Fig. [T] 
middle and Fig. |2l bottom). With its flux density of ~ 90 mJy, 
it was clearly the brightest radio source in the field (at least 
five times brighter than the second brightest source). Almost 
one year later, on 2009 April 08, the source had faded to a 
flux density of ~ 11 mJy, i.e. losing almost 90% of its flux 
density (see Fig. [T] bottom). We estimated its position to be 
Q'y20()o=09''55'"5 1 .55 1 ' + 0.008' 5y20oo=69°40'45.792" ± 0.005". 
Following the traditional convention of previous papers nam- 
ing a source after the position offset to a;Bi95o=09''55"' and 
(Jgi9g,)=69°40', our new detection would be 42.82+59.54. The 
positions in the two observations are consistent and the un- 
certainty was estimated by comparing the position of another 
clearly identified source (44.0+59.6) with the positions from ob- 
servations with the Mul ti-Element Radio Lin ked Interferometer 
Network (MERLIN) bv lMuxlow et"!!] ( |1994|) . 

To constrain the start time of the flare, we reduced two earlier 
VLA observations of M82. The spectral line observation on 2008 
March 24 had lower quality, since we observed with a narrow 
bandwidth and had only ~6 minutes integration time on M82. 
Nevertheless, we could easily confirm the detection of the new 
source at the same position with a similar but slightly higher 
flux density as in the 2008 May 03 observation. Next, we ana- 
lyzed 4.8 GHz data from 2007 October 29. Here, no bright point 
source was discovered (see Fig.[T] top), so we conclude that the 
start of the flare was between 2007 October 29 and 2008 March 
24. However, it is possible that the source could be highly self- 



absorbed early in its development and thus might have been de- 
tectable at high frequencies, but not at 4.8 GHz. In this case, the 
flare might have started earlier, but not earlier than 2007 January 
28 (the day of our first 22 GHz observation, see Fig.|2l top). 

The comparision of flux densities from different epochs in 
M82 is difficult due to the diffuse emission, in particular when 
comparing observations with different spatial resolutions. In 
Table [1] we list the flux densities of the new source along with 
the flux densities of one additional source in M82 (44.0+59.6), 
M81*, and three background quasars. 

4. Pre-f lare observations at other wavelengths 

iMatsumoto et alj (1200 Ih present high-resolution (FWHM ~ 
0!'5) x-ray imaging of the central 1' x 1' (1.1 x 1.1 kpc) re- 
gion of M 82 with the High-Resolution Camera (HRC) aboard 
the Chandra X-ray Observatory. These images, taken on 1999 
October 28 and 2000 January 20, show a total of 9 sources 
within this area, some of which are highly variable. They 
do not find a counte rpart to our new radio source. Neither 
do iKong et al] (|2007|) in a total of 12 datasets of a similar 
sized region taken with the Chandra HRC and Advanced CCD 
Imaging Spectrometer Array (ACIS-1 and -2) taken between 
1999 September 20 and 2005 August 18. These authors also 
present Hubble Space Telescope (HST) //-band (1.6 fj.m) imag- 
ing of the region with the Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) that, again, does not show a 
co unterpart to our variable radio source. 

iKording et all (|2005) observed M82 with the VLA at 8.4 
GHz eight times between 2003 June and October in order to 
detect possible radio flares from several ultraluminous X-ray 
sources (ULX). No significant emission was found at the po- 
sition of our new detection above a noise level of 70 juJy at each 
epoch. A very sensitive eight-day integration of M82 at 5 GHz 
wavelengt h was performed w ith MERLIN between 2002 April 
1 and 28 jFenech etalj|2008 ). The resulting 40 millarcsecond 
resolution images show no significant emission at the position 
of our new sourc e above the rms noise level of 17 //Jy beam 
iTsai et alj (l2009l) present a 7 mm map from the VLA from 2005 
April 22, with no detection of more than 0.3 mJy at the position 
of our source. 

5. Discussion 

The most straightforward explanation for this new source is a 
new radio supernova. A quantitative analysis of the lightcurve is 
difficult due to the small number of data points (3) we obtained 
with different angular resolutions, the complication of diffuse 
emission in M82, and the uncertain absolute flux scale in our 
observation of 2008 May 03. We extracted the flux d ensities of 
our source by restricting the interferometer (u,v)-data to >30 kA 
to remove most of the extended emission (Table[T]). 

We fitted the lightcurve of the source with an exponen- 
tial decay: S(t)=S() e*'""''^^''. This yields a decay timescale of 
T,/=0.46+0.03 yr (Fig.[3]l. Since we do not know the exact time 
of the onset of the flare, to, we get values of So between 1 10 and 
270 mJy for to between 2008 March 24 and 2007 October 29. 
For the fit, we added in quadrature the difference between peak 
and integrated flux density as one error estimate and systematic 
(flux density scale) errors of 5% for the first and third epochs, 
and 10% for the second epoch, where no proper flux density 
calibrator was observed. The resulting was 0.5. Based on 
this fit, we predict that the source flux density will drop to 5 
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Table 1. Details of the detected compact sources in our VLA ob- 
servations with statistical errors on the 22.2 GHz flux densities. 



Source 


Date 


Peak Flux 


Integrated Flux 






[mjy beam"'] 


[mJy] 


M82 transient 


24/03/2008 


104.4±1.3 


99.6+2.3 


(42.82+59.54) 


03/05/2008 


89.9+0.1 


88.4+0.2 




08/04/2009 


11.1+0.1 


9.2+0.2 


44.0+59.6 


28/01/2007 


15.0+0.2 


77.8+1.0 




03/05/2008 


12.5+0.1 


11.0+0.2 




08/04/2009 


14.7+0.1 


12.4+0.2 


M81* 


28/01/2007 


71 + 1 


72+1 




03/05/2008 


150 


150 




08/04/2009 


195+3 


179+5 


0945+6924 


28/01/2007 


14.3+0.2 


14.4+0.3 




03/05/2008 


14.8+0.1 


14.7+0.2 




08/04/2009 


13.8+0.2 


12.2+0.3 


0948+6848 


28/01/2007 


22.5+0.2 


24.3+0.3 




03/05/2008 


59.1+0.2 


59.8+0.3 




08/04/2009 


63.6+0.4 


56.3+0.6 


1004+6936 


28/01/2007 


31.0+0.2 


30.8+0.3 




03/05/2008 


30.9+0.2 


30.4+0.2 




08/04/2009 


34.6+0.3 


30.8+0.4 



mJy by mid 2009 and 1.5 mJy in early 2010. We also fitted a 
simple exponential decay to the 22 GHz lightcurve of SN199J 
in M8 1 . Here we used the first year of data a fter the peak in 
the lightcurve published in IWeiler et al.l ( |2007|) . The fit yields 
T^=0. 72+0.09 yr, indicating that the new source in M82 decays 
faster than SN1993J. 

A power-law fit (S(t)oc(to-t)") is not consistent with the 
lightcurve of the new source in M82 Ct^^y = 5-8, for to be- 
tween 2007.8 and 2008). If one allows an earlier to (e.g., if the 
source is highly self-absorbed at lower frequencies), the power- 
law fit improves (x^p^f - 2, for to=2007.08, one day after our 
2007 January 28 observation). 

The position of the new source is marginall y consistent with 
the p osition of the kinematic center of M82 (IWehachew et alj 
119841) . This raises the possibility that, rather than a supernova, 
we could have detected a flare from a supermassive black hole 
in the center of M82. Flar es in AGN sources c an often be fitted 
with exponential decays jValtaoia et alJlT999h . For example, a 
strong radio flare in 1999 in the Seyfert galaxy II IZw2 has a 
decay rate of tj=0.73 year (iBrunthaler et al.ll2005ah . 

Other explanations, such as luminous flares from quies- 
cent supermassive black holes induced by a close passage of 
a star that is torn apart by tidal forces, are also possible. 
However, such flares show a power-law deca y with a - -5/ 3 
d^ans & Kochanek 1989; Ayal et al. 2000; Ge zari et al.|[2009l) . 
which is not consistent with our lightcurve. 

However, M82 has never shown evidence of a nuclear su- 
permassive black hole (which would be surprising for a small 
irregular galaxy). Since the progenitor for our flare showed no 
X-ray emission, a stellar or intermediate mass black hole is also 
not probable. Thus, based on the current data, a new radio super- 
nova seems to be the most likely explanation. 
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RIGHT ASCENSION M2000t 



Fig. 1. Top: VLA B-configuration image at 4.8 GHz of M82 
on 2007 October 29. Contours start at 3 mJy and increase with 
factors of 2. The beamsize is 1.7x1.1 arcseconds with a position 
angle of -58°. There is no detectable source at the position of 
our new source. The peak flux in the image is 45 mJy beam"'. 
Middle: VLA C-configuration image at 22 GHz of M82 on 2008 
May 03. Contours start at 0.7 mJy and increase with factors of 
2. The beamsize is 0.90x0.83 arcseconds at a position angle of - 
73°. The new source at position a'y20oo=09''55"'51.55r ± 0.008* 
572000=69°40'45.792" ± 0.005" is cleai-ly visible, with a peak 
flux of 90 mJy beam ' . Bottom: VLA B-configuration image at 
22 GHz of M82 on 2009 April 08. Contours start at 0.7 mJy and 
increase with factors of 2. The beamsize is 0.32x0.29 arcseconds 
with a position angle of -60°. The new source is still clearly vis- 
ible, but its peak flux density has decreased to 1 1 mJy beam"'. 
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Fig. 2. Contours and grey scale represent 22.2 GHz images of the 
central region of M82 taken on 2007 January 28 (top) and 2008 
May 03 (bottom). The new bright source is very conspicuous 
in the latter image. Contour values are -4, 4, 8, 16, 32, 64, 128, 
256, and 512 times 0.2 mJy beam"', roughly the (comparable) 
rms noise level in both images. The sta r symbols mark the posi- 
tions of the x-ray sources discussed bv lMatsumoto et al] (1200 ll) 
and in the top panel bear these authors' nomenclature. The as- 
terisk gives the position of the radio kinem atic center (R.K.C.) 
of M82 as determined bv lWeliachew et akl (il98 4). Both images 
were restored with a circular beam of 1.5 arcseconds. Since the 
observation on 2007 January 28 was made in C configuration, 
there was much more data from short baselines, which results in 
better sensitivity to extended emission. 
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Fig. 3. Lightcurve of the new radio transient from our observa- 
tions (squares). Also shown are an exponential-decay fit (solid 
Une) and power-law fits (dashed fines), assuming difi'erent peak 
times of the flare: 2008.0 (lower line), 2007.8 (middle line), and 
2007.08 (upper fine). 



